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term.  Specifically  for  plane  wave  incidence  on  a  Gaussian  distributed  rough  surface 
with  a  Gaussian  correlation  function  we  anaiyticadly  evaluate  the  incoherent  intensity 
for  single  scattering  with  correlated  heights  and  slope,  and  for  double  scattering  where 
the  heights  and  slopes  are  uncorrelated. 
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ABSTRACT 


^  Rigorous  analytical  methods  in  rough  surface  scattering  lead  to  integral  equations. 
In  coordinate  space  these  integral  equations  are  written  on  the  surface  value  of  the  field, 
(Neumann,  hard,  TM  boimdary  value  problem), ''the  surface  value  of  the  normal 
derivative  of  the  field  (Dirichlet,  Soft,  TE)^pr  coupled  equations  on  both  (interface  or 
penetrable  surface  or  dielectric  problem).  An  additional  integration  is  necessary  to  find 
the  scattered  or  transmitted  field.  Li  Fourier  transform  space  the  corresponding 
equations  are  on  the  scattering  and/or  transmission  amplitudes  with  a  direct 
interpretation  as  scattered  or  transmitted  fields  as  a  function  of  angle.  -  i  "  a  .:  j- 

The  integral  equations  express  the  full  multiple  scattering  solution  and  this  must 
be  treated  numerically.  Truncation  of  the  equations  can  lead  to  analytical  results  for 
single  scattering  (which  is  well  known)  and  for  double  scattering.  We  describe  here  the 
approximations  necessary  to  analytically  evaluate  to  closed  form  the  double  scattering 
term.  Specifically  for  plane  wave  incidence  on  a  Gaussian  distributed  rough  surface 
with  a  Gaussian  correlation  function  we  anadytically  evaluate  the  incoherent  intensity 
for  single  scattering  with  correlated  heights  and  slopes,  and  for  double  scattering  where 
the  heights  and  slopes  are  imcorrelated. 
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1.  INTRODUCTION 


Several  years  ago  we  introduced  Jb-space  integral  equations  for  scattering  from  a 
rough  surface  using  Feynman-diagram  and  partial  summation  techniques  [1-3]. 
Numerical  solution  of  a  truncated  integral  equation  for  the  scattering  amplitude  of  the 
coherent  wave  using  Gaussian  statistics  on  the  surface  illustrated  the  necessity  of  using 
multiple  scattering  techniques  of  this  type  to  describe  the  scattering  [4].  For  the 
coherent  field  the  results  agreed  well  with  experimental  data  into  the  regime  of  large 
roughness.  A  review  of  the  specific  work  as  well  as  a  review  of  the  field  can  be  found  in 
DeSanto  and  Brown  [5]. 

One  approximation  used  in  this  previous  work  was  related  to  the  specific  way  the 
transform  equation  arose  from  the  coordinate-space  equation.  Recent  work  on 
backscatter  enhancement  [6,7]  has  prompted  us  to  attempt  to  remove  this 
approximation  in  order  to  have  an  exact  version  of  the  ^space  equation,  and  to 
thereby  enable  us  to  study  the  single  and  double  scattering  approximations  of  the  exact 
result.  The  latter  define  the  boundary  between  possible  zuialytical  approximations  and 
necessary  numerical  results  in  the  sense  that  triple  and  higher  order  scattering  must  be 
done  numerically.  We  briefly  describe  the  successful  result  of  this  attempt  for  the 
Neumann  or  hard  boundary  value  problem  here. 

2.  PHASE  MODULATION 

For  a  plane  wave  incident  at  angle  from  the  normal,  the  phase  change  (with  no 
horizontal  translation)  upon  reflection  into  angle  B,  from  the  rough  surface  z  =  h{xf) 
(with  zt  =  (i,  y))  is  given  by 


—  hkcosBi  —  hkcosBf  =  {kg  — kg)h  ,  (1) 

where  k  is  the  wavenumber  with  incident  z-component  k\  =  —k  cos  0,-  and  scattered  z- 
component  kg  =  kcoaBf.  With  an  additional  horizontzd  translation  the  full  wave 
scattering  fvmction  can  be  written  as 


exp 


(2) 


where  k‘(ifc)  is  the  incident  (scattered)  direction  and  Z/,  =  (z,  y,  h)  is  a  point  on  the 
siirface.  If  we  multiply  this  by  the  (non-unit)  normal 


«m(a^)  =  ^m3  -  h{xt)  ,  (3) 

where  is  the  transverse  difierential  operator,  and  Integrate  over  the  trzuisverse 
coordinates  we  get  the  fimction 
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Am{k-k') 


(4) 


which  is  termed  the  phase  modulated  amplitude  spectrum.  It  describes  all  the 
scattering  from  initial  state  k*  into  scattered  state  k.  It  is  also  just  the  physical  optics 
approximation  for  the  scattering.  It  will  occur  again  below  in  describing  the  multiple 
scattering  in  the  integral  equation  and  of  course  in  the  single  and  double  scattering 
temos. 


S.  GREEN’S  THEOREM  AND  THE  SCATTERING  AMPLITUDE 

Using  Green’s  theorem  the  field  scattered  from  a  rough  surface,  ip*,  can  be  written 
at  a  point  x  as 

=  // G{x,XnW{xh)  -  G(i,x;)jv*(i;)|dz;  (s) 

in  terms  of  its  value  on  the  surface  and  the  value  of  its  normal  derivative  on  the  surface 


N‘{xh)  =  nmdm^*{xh) 


(6) 


where  repeated  subscripts  are  summed  from  one  to  three.  Here  G  is  the  free  space 
Green’s  function 


G{?,  x') 


exp 


iko  (  I  —  i' 


47r  1 1  —  i' 


(7) 


md  n'„  is  just  n„{xt).  The  asymptotic  value  of  the  Green’s  function 
(J?  =  \x\,x  large) 


.UfoR  , 


yields,  using  (5),  the  asymptotic  value  of  the  scattered  field 

ik^R 

n.)  ~  ’ 


(8) 


(9) 


where  F{k)  is  the  scattering  amplitude 
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^(^)=-//«  ’*  **  dx'i  .  (10) 

It  appears  as  a  quasi-Fourier  transform  over  the  surface  field  and  normal  derivative. 
The  object  is  to  find  an  integral  equation  on  F  and  determine  approximations  from  it. 
The  actual  scattering  amplitude  is  found  from  (10)  set  on  the  energy  shell,  i.e.  set 

K  =  (*o  • 

4.  NEUMANN  PROBLEM 

The  Neiimann  (N)  boundary  value  problem  is  given  by  the  boundary  condition 

=  n^a„,^^(i/,) 

=  -nm^m^'(?A) 

=  -  N\xh)  (11) 

which  states  that  the  normal  derivative  of  the  total  field  (V'*  +  V>*)  vanishes  on  the 
surface.  Here  V'*  is  the  incident  field.  The  result  can  be  written  from  (10)  as 

F''W  =  F^(k)-iK,K(!s)  .  (12) 

where 

^f(*)  =  f  f  ,  (13) 

and 

0m(*)  =  //«’***  ,  (14) 

where  will  be  shown  to  describe  single  scattering  (and  is  known  since  the  incident 
field  is  known)  and  where  the  object  will  be  to  find  an  integral  equation  on  th® 

remaining  term  in  (12). 

In  order  to  do  the  latter  it  is  first  necessary  to  formulate  an  integral  equation  on 
0*.  The  latter  follows  from  the  limit  of  (5)  as  i.e.  as  x  approaches  a  point  on 

the  surface,  and  the  use  of  (11).  The  normal  derivative  of  G  in  (5)  has  a  discontinuity 
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on  the  surface,  and  the  regularized  integral  equation  can  be  written  as 


=2// x'k)N'{xh)dxt  -  2  J  j n‘„  Rm{xh  -  Xh)rl)*{xh)dx't  ,  (15) 


which  is  a  second  kind  equation  with  a  regularized  kernel  gi^en  by  (G  is  the 
Fourier  transform  of  G) 


Rm(j)  =  M-’/ //«*■•  G(k)P„{k)dk  ,  (16) 


where 


Pm{k)  =  t 


^  +  ^m3 


(17) 


with  —  kf.  The  symbol  P  represents  the  Cauchy  principal  value  distribution. 

We  can  then  “Fourier”  transform  (15)  using  (14)  to  find  am  equation  on  (i>m  and 
thus  satisfy  (12).  The  result  can  be  written  in  the  form 

Fnik)  =  f!' (k)  +  F? {k)  +  {k)  ,  (18) 


where  Fj  is  given  by  (13),  and  F^  (k)  is 

F?(k)  =  -2ik„ / / c-*  ‘‘n„[StP){z',)dz;  ,  (19) 

where  S  N*  vs  the  single  layer  potential  with  density  N*  given  by 


{SIir)(z)  =  J  J  G{z,zi,)  N‘{zS)dzi  . 


(20) 


Given  the  incident  field  this  function  is  also  known.  Its  value  on  the  surface  is  just  the 
continuous  limit  of  its  value  as  the  field  point  x  approaches  the  surface  (x-+ 1/,). 

Finally  (k)  is 


F^l.k)  =  2ik„\2z]-^J  J  J  A^(k-g)G{fl)  W(g)de  ,  (21) 
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where  W  satisHes  the  integral  equation 


i*'(s)=  »''(a -2(2irr’p„  (^///-  m{S-Q')G{nn&')dQ  ,  (22) 

whose  Bom  term  is 

W‘  iS)  =  2PM  J / c"‘'^'*‘n;„(5iV){x;)dx;  ,  (23) 

which  equals  i  {Q)  on  the  energy  shell  =  (^  ~  j  • 

5.  PLANE  WAVES 

For  plane  wave  incidence 


V'*(f)  =  exp(ii* 'x)  , 


(24) 


we  have,  using  (13),  that 


.  (25) 

which  describes  single  scattering  (one  A-function)  from  incident  state  k*  to  scattered 
state  k.  F2  can  be  written  using  (19)  as 


F?(k)  =  2\2^]-^k„ki  J  J  jA„(k-g)  G(0)A,{g-k‘)dS  ,  (26) 


which  describes  double  scattering  (two  A-functions)  from  incident  state  k*  into  all 
states  propagation  in  these  states  (  with  the  propagator  G  ),  and  rescattering  into 
state  k.  This  is  just  the  matrix  element  of  the  double  scattering  term  between  an 
incident  plane  wave  state  and  scattered  state  k.  The  function  F^  in  lowest  order  has 
three  A-functions  and  describes  triple  and  higher  order  scattering.  The  functions  Fi 
and  F2  appear  to  be  the  limiting  case  for  analytical  approximation  methods.  Triple 
and  higher  order  scattering  must  be  approached  numerically.  Note  that  in  terms  of  the 
separation  we  described,  both  single  and  double  scattering  terms  are  given  explicitly  in 
terms  of  the  incident  field.  This  is  also,  of  course,  true  for  the  more  general  case  of 
non-plane  wave  incidence,  zdthough  for  the  latter  more  integrals  are  involved. 
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6.  KIRCHHOFF  APPROXIMATION 

It  can  easily  be  seen  that  the  (single  scattering)  Kirchhoff  approximation  for  plane 
wave  incidence  applied  to  (10)  yields  the  result 

=  -<.»(*-*■)  •  (27) 

Comparison  with  the  exact  single  scattering  result  (25)  indicates  that  the  Kirchhoff 
approximation  gets  some  of  the  single  scattering  right,  but  then  approximates  the 
remaining  true  multiple  scattering  by  a  single  scattering  term  involving  the  scattering 
state  ib. 

7.  AVERAGES 

There  are  two  ways  to  determine  average  field  or  intensity  quantities  if  the  surface 
is  random.  One  is  to  do  a  straightforward  numerical  solution  of  the  deterministic 
equations  for  each  member  of  the  statisticzil  ensemble  of  surfaces  (Monte  Carlo 
calculation),  and  then  to  form  the  averages  of  the  desired  quantities.  For  a  one¬ 
dimensional  siirface  this  has  been  done  in  coordinate-space  by  Thorsos  [8],  The  method 
includes  all  orders  of  multiple  scattering  since  it  is  the  solution  of  an  integral  equation. 
The  second  way  is  to  apply  averaging  operators  to  the  equations  to  generate  integral 
equations  on  the  field  average  (coherent  field  or  Dyson  equation)  or  the  average  of  the 
second  (or  higher  order)  moments  of  the  field  (Bethe-Salpeter,  etc.  equations).  In  k- 
space  and  using  diagram  methods  we  have  previously  described  this  procedure  [1-3,5]. 
Here  we  adopt  the  latter  philosophy  but  only  applied  to  explicit  representations  of 
single  and  double  scattering. 

The  averaging  operator  E  k  an  integral  operator  whose  dimensionality  (and 
integration  variables)  is  the  number  of  random  variables  involved  and  whose  kernel  is 
the  joint  probability  density  fimction  of  all  the  random  variables.  We  have  expressed 
the  scattering  amplitudes  in  terms  of  A,„-functions,  and  each  A„,-function  contains 
three  random  variables,  the  height  function  h,  and  the  two  slope  terms,  and  hy. 
Coherent  averages  are  defined  (here)  as  first  taking  the  averages  of  the  scattering 
amplitudes  Fi  or  F2  and  then  forming  the  coherent  intensity  from  the  square  of  the 
magnitude  of  this  result.  For  coherent  single  scattering  we  thus  have  three  random 
variables,  and  for  coherent  double  scattering,  six.  For  Gaussian  height  distribution  and 
Gaussian  correlation  fimction  (or,  in  general,  for  a^iy  statistically  homogeneous  or 
spatially  stationary  process)  it  can  be  shown  that  the  coherent  averages  for  any  order 
of  multiple  scattering  yield  intensities  which  are  non-zero  in  only  the  specular  direction. 
For  incoherent  scattering  (here)  we  first  form  the  magnitude  squaired  of  the  scattering 
amplitudes  and  then  apply  the  averaging  operator.  For  incoherent  single  scattering  we 
thus  have  six  random  variables  (two  A-functions)  and  for  incoherent  double  scattering, 
twelve.  The  latter,  of  course,  begins  to  get  formidable. 
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In  spite  of  the  apparent  multivariability,  the  restriction  to  Gaussian  statistics 
results  in  only  two-point  moments  (correlation  function)  on  the  surface  heights  and 
slopes.  There  is  only  one  correlation  function  but  because  of  the  multivariability  it  has 
several  combinations  of  coordinate  variability.  Many  coordinate  derivatives  also  occur 
(for  height-slope  correlations)  and  second  derivatives  (slope-slope  correlations).  It  is 
tempting  becaiise  of  aJl  this  variability  to  cany  out  a  partial  integration  of  the  A^- 
functions.  The  result,  neglecting  the  partially  integrated  term  is 

a  k„(k,)-‘ A[k)  ,  (28) 


where 


(29) 


and  thus  reducing  the  random  variability  of  each  A^-fnnction  to  a  single  random 
variable,  the  surface  height.  This  essentially  deconelates  heights  and  slop>es. 
Attractive  as  this  approximation  is  in  terms  of  analysis,  there  is  no  tmderstanding  of  its 
effect  in  terms  of  physical  modelling,  in  particular  for  scattering  processes  involving 
large  heights  and  slopes. 

8.  EXAMPLES 

For  a  Gaussian  process  we  list  two  examples  of  the  incoherent  intensity.  One  is 
for  single  scattering  and  one  for  double  scattering.  For  single  scattering  the  incoherent 
(7)  intensity  is 


(30) 


where  for  we  give  the  result  in  such  a  way  that  it  includes  the  exact  single 
scattering  result  (25)  or  the  Kirchhoff  result  (27).  It  is  a  single  integral  and  includes  adl 
height  slope  correlations 


00  5 
^i;(*)  =  J  pexp  -  (koca  (1-m(p))]  ^Vj{0,p)dp  , 


(31) 


where  Is  the  ensonified  surface  area,  L  Is  the  correlation  distance,  ko  the 
wavenumber,  a  the  rms  height,  c  =cos^ +  cos^*,  pi{p)  the  Gaussian  correlation  function 
(  where  p  is  dimensionless)  and 
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(32) 


VjCi?,  p)  =  Apn[p)Ji{Mp)  ckQa{(jl  L){KxCosB  + K2smB)K:i  , 


V:,{e,p)  =  2n{p)Jo{Mp){a/L)\K\+Kl)  , 


^4(^>p)  = -4p*m(p)  l  +  (fcoca)*/i(p)j  Jo(A^^p)  ^  ifi  cos^  +  ^^2  sintf  (35) 


P)  =  4M  ^p/i(p)  1  +  (A;o  ca)^/i(p)j  J i  (Mp) 


(a/I)  2j(iiricos(9  +  iir2sin6>)2  -  jiKl  +/ri)i  , 


where  the  Ts  are  Bessel  functions  and 

M=\kt~k\\L  .  (37) 

Here  if  Kj  =  k'j  —  kj  we  have  the  Kirchhoff  result  from  (27)  and  if  Kj  =  fcj-  the  exact 
single  scattering  result  from  (25).  With  a  Taylor  expansion  in  the  exponential  the 
integrals  can  be  evaluated  using  tables. 

The  second  example  is  that  of  exact  double  scattering  for  a  Gaussian  process 
where  in  order  to  carry  out  the  averaging  it  was  necessary  to  cairry  out  a  partial 
integration  of  the  A^-fimctions  using  (28)  (  i.e.  heights  and  slopes  are  not  correlated). 
In  addition  the  Taylor  expansion  mentioned  to  evaluate  (31)  is  used.  The  result  is 


where 


A;= 
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A|  =  W-W/tti-y  .  (40) 

(which  is  a  transverse  vector) , 

^•'  =  2(l  +  A?)-‘(*'A<-ti-A,)  (42) 

<•  =  (kif  (1  +  A?)-  ‘  (A?  - 1  -  2  tj  •  A,/*  ')  (43) 


and  S  is  just  5*  with  /b*  replaced  by  k.  Although  it  has  been  possible  to  do  the  full 

evaulation  it  is  under  the  restriction  of  uncorrelated  heights  and  slopes.  I  regard  this  as 

a  severe  restriction  and  am  presently  attempting  to  overcome  it. 
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